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Abstract 

The equation of state (EOS) of finite nuclei is constructed in the rela- 
tivistic Thomas-Fermi theory using the non-linear a — lo — p model. The 
caloric curves are calculated by confining the nuclei in the freeze-out volume 
taken to be a sphere of size about 4 to 8 times the normal nuclear volume. 
The results obtained from the relativistic theory are not significantly differ- 
ent from those obtained earlier in a non-relativistic framework. The nature 
of the EOS and the peaked structure of the specific heat C v obtained from 
the caloric curves show clear signals of a liquid-gas phase transition in finite 
nuclei. The temperature evolution of the Gibbs potential and the entropy 
at constant pressure indicate that the characteristics of the transition are 
not too different from the first-order one. 
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1 Introduction 



Nuclear caloric curves have been obtained from a number of experiments in recent 
times from energetic nucleus-nucleus collisions. In the GSI data [|lj for Au + 
Au at 600 AMeV, the temperature was found to be practically constant at a 
value of T ~ 5 MeV over the excitation energy (e*) range of 3 - 10 MeV per 
nucleon after which e* was found to increase linearly with temperature as in a 
classical gas. This is suggestive of a sharp liquid-gas type phase transition. The 
caloric curve obtained from the collision of Au on C at 1 AGeV in the EOS 
collaboration experiment also shows a plateau at T ~ 6 MeV; this is not as 
prominent as in the earlier case. However, the heat capacity C v derived from 
this caloric curve of the EOS group shows a peaked structure at T ~ 6 MeV 
indicating existence of a phase transition. Even at a relatively low bombarding 
energy of 47 AMeV j3| for several reactions, it has been seen recently that the 
caloric curves show plateau at T ~ 7 MeV in the excitation energy range of ~ 
3.5 to 7 MeV per nucleon. It would be interesting to know whether this is a 
precursor to the liquid-gas phase transition in the finite nuclei. Theoretical analysis 
of infinite nuclear matter (symmetric as well as asymmetric), both in the non- 
relativistic 0, [5|] and relativistic framework || [71 predict Van der Waals type 
isotherms in their equation of state (EOS) implying coexistence of liquid and gas 
phases. Finite size effects and the Coulomb interaction between protons might 
change such a behavior, but the EOS of realistic nuclei, calculated only recently || 
in the non-relativistic Thomas-Fermi (TF) theory display same kind of isotherms, 
liquid-gas coexistence and a liquid-gas phase transition at a temperature quite 
below the critical temperature for infinite nuclear matter. The calculated transition 
temperatures are still somewhat higher than the observed ones, the caloric curves 
do not match exactly the ones derived from experiments, but the calculations in 
Ref . give unmistakable signals of the liquid-gas phase transition in the finite size 
nuclear systems. 

Relativistic mean field (RMF) theories have been applied successfully to explain 



the ground state properties of nuclei over the entire periodic table [ TO | . This theory 
has also proved to be very fruitful in explaining various details of exotic nuclei near 



the drip lines pT| , |12[ . In contrast to the non-relativistic models, the RMF theory 
uses a single set of parameters to explain all these properties. It would therefore 
be very interesting to investigate the EOS of finite nuclei and the related exotic 
phenomena like liquid-gas phase transition in the relativistic approach. The present 
paper aims at understanding these thermodynamic properties of hot nuclei in a 
relativistic Thomas-Fermi (RTF) theory. 

In Sec. 2, we briefly outline the formalism used. The results and discussions 
are presented in Sec. 3. The summary and conclusions are given in Sec. 4. 



2 Formalism 

A brief outline of the calculations of the relevant thermodynamic quantities in 
the relativistic Thomas Fermi approximation is presented in this section. The 
Lagrangian density used is given by 



£ = %(i^-M)% + ^(7d tl a-U((7)-g a %(7% 

-gjbcfp^i - \f^f^ - e^ (1 ~ T3) A^. (1) 

The meson fields included are those of the isoscalar a meson, the isoscalar- vector uj 
meson and the isovector- vector p meson. The arrows in Eq. (|XJ) denote the isovector 
quantities. The z-component of isospin, r 3 , is taken to be +1 for neutrons and —1 



for protons. For appropriate description of the nuclear surface properties ||13 
scalar self-interaction term U(a) of the a meson is included in the Lagrangian, 

U(a) = l -mla 2 + ^ 2 a 3 + -^a\ (2) 

The quantities M, m a , m w and m p are the nucleon, a, uj and the p— meson masses, 
respectively, while g a , g^, g p and e 2 /47r = 1/137 are the corresponding coupling 



constants for the mesons and the photon. The field tensors of the vector mesons 
and of the electromagnetic fields have the following structure: 

W = (3) 

$r = &*F-drF-g p (pxF), (4) 

f" = d» A v - d u A» . (5) 

The equations of motion are obtained from the variational principle. The mean 
field approximation is introduced at this stage by treating the fields as c— numbers 
or classical fields. This results in a set of coupled equations, namely the Dirac 
equation with potential terms for the nucleons and the Klein-Gordon type equa- 
tions with sources for the mesons and the photon. Time reversal invariance and 
charge conservation get the equations simplified in the static case. The resulting 
equations, known as relativistic mean-field equations, have the following form. The 
Dirac equation for the nucleon is 



where 



and 



{-ia • V + V(r 3 , r) + f3 [M + <S(r)]}^ = e^, (6) 



(1 — T ) 

V(r 3 , r) = 9u u (t) + g p T 3 p (r) + e l - 3) A (r), (7) 



5(r)=^(r), (8) 

are the vector and the scalar potentials, respectively. The scalar potential con- 
tributes to the effective mass as 

M*(r) = M + S(r). (9) 

The Klein-Gordon equations for the mesons and the electromagnetic fields with 
the nucleon densities as sources are 

{-A + ml}a(r) = -g aPs (r) - g 2 a 2 (r) - g 3 a 3 (r), (10) 

{-A + ml}u (r) = g uPv (r), (11) 

{-A + mj}p (r) = g p p 3 (r), (12) 

-AAo(r) = ep c (r). (13) 



While considering a finite nucleus, for simplicity, we assume it to be spherically 
symmetric. The above field equations then can be written in a general form 

rrirj) are the meson masses for <ft = <j,lj and p and is zero for the photon. The source 
term S^r) is given by the right hand side of the Eqs. ( |T0| ) - C p~3| ) for a,uj,p and 



Coulomb fields. The above equation (|HD can be solved using Green's function [ |T0f 

0(r) = / r^^r^S^dr', (15) 



where 



G^r, r') = L_[ e -™^-V| _ g-m^lr+r'^ 

for the massive fields and 

{- for r > r' 

; (17) 
4 for r < r', 

for the Coulomb field. 

The quantities p s , p v , p% and p c appearing on the right hand side of Eqs. ( [TOp - 



( |13D are the scalar, baryon, isovector and charge densities, respectively. They can 



be obtained as 



Ps{r) = £p s (r 3 ,r), 

T3 



i ll 

T3 

Pc(r) = £(^)p.(T 3 ,r). (18) 

In the Thomas-Fermi approximation, the quantities p v {^3, r) and p s (T3, r) are given 
as 

7 f°° *t, 



ft,(r 3 ,r) = / f(e,T)k 2 dk, (19) 

Z7T 







PsM = ZL= f(e,T)k?dk, (20) 

2tt 2 Jo Jtf + M* 2 (r) 



where the spin degenaracy factor 7 is equal to 2. The self-consistent occupancy 
function /(e, T) is obtained through the minimisation of the thermodynamic po- 
tential 



G = E — TS — fiN, (21) 

and is given by 

^■ r )= 1 + e L)/T . M 

with 

e(r 3 , fc, r) = V(r 3 , r) + + M* 2 (r). (23) 

The chemical potential \i is adjusted to get the desired number of particles (neu- 
trons and protons) given by 

n(r 3 )=A7i r 2 p v {T 3 ,r)dr. (24) 
Jo 

In Eq.(24), r max determines the confining volume V taken to be spherical. From 
Eqs.(19) and (22), it can be seen that at large distances the baryon density p v ~ 
e (p-M)/ T a nd is therefore a non-zero constant at finite temperature. The solution to 
baryon density and hence the various observables depend on the choice of the size 
of the box in which the calculation is performed. At zero temperature, however, 
the solution is independent of the choice of the confining volume once it is larger 
than the normal nuclear volume Vq. Exactly the same characteristic is seen in the 
nonrelativistic case [14, [0|. The choice of the volume for the evaluation of the 



thermodynamic variables would be discussed in the next section. 

For a nuclear system with mass number A, the total energy E(T) is given by 

0, 

E(T) = E part + E a + E uNL + E u) + E p + E c + E CM - AM, (25) 

with 

2 7 



V^y rTmax POO 

E part = — E/ rdr k 2 e(r 3 ,k,r)f(e,T)dk, (26) 

7T ~ JO JO 

E a = ~g a J d 3 rp s (r)a(r), (27) 



EaNL — 


1 r n ( 1 o 1 1 

-i/A{i^(r) + I^(r)}, 


(28) 


E„ = 


~2 9 " J rf3r A>( r V°( r )' 


(29) 


E P = 


~2^/ d 3 rp 3 (r)p°(r), 


(30) 


E c = 


/ d'rp c (r)A°( r ), 


(31) 


Ecm = 




(32) 



The free energy F is given by (E — TS 1 ) where the entropy S can be calculated 
from the Landau quasiparticle approximation, 

S = r 2 dr k 2 iflnf + (l-f)ln(l-f)}dk. (33) 

7T " JO JO 

The specific heat , C v , and the pressure, P, can be calculated from 

dE 



Cv <rr 



p = 

dv 



(34) 
(35) 

T 



The baryonic density, the mesonic and the Coulomb fields are obtained iteratively 
through the following scheme: 

i) An initial guess is made for the fields c(r), ou(r), p(r) and A (r). 

ii) The effective mass M* and the energy e given by Eqs. (9) and (23), respectively 
are calculated with the guess fields. The proton and the neutron chemical poten- 
tials (p,) are adjusted to reproduce the given number of nucleons of each kind. 

iii) From Eq.(18), the various densities and hence the source terms are calculated. 

iv) These source terms are used in Eqs. (10)- (13) for the generation of the new 
fields. 

The steps (ii)- (iv) are repeated till the desired accuracy is reached. 

The expressions for the EOS for nuclear matter can be obtained as a special 
case of a finite nucleus by ignoring the gradient terms in the field equations. This 
simplifies the expressions for the relevant observables and are given in Ref . . 



3 Results and discussions 



In this section, the results of our calculations for the EOS of infinite symmetric 
and asymmetric nuclear matter as well as of a few finite nuclear systems are first 
given. We have chosen 40 Ca, 109 Ag and 150 Sm as the reprersentative systems. The 
results of our calculations of the caloric curve for these nuclei are next presented. 
A host of parameter sets for the nonlinear o — uj — p model are available which 
produce almost similar ground state properties of nuclei over the whole periodic 
table but with widely different values of nuclear incompressibility (Koo). To study 
the effects of different on the results of our calculations, we have chosen the 
parameter sets NL1, NL3 and NLSH [II], having equal to 212, 272 and 356 
MeV, respectively. 



3.1 The nuclear EOS 

The EOS of symmetric and asymmetric nuclear matter in the RTF theory has al- 
ready been calculated by Mwller and Serot |7j . For comparison of results for finite 
nuclei obtained with a given set of parameters and also for completeness, we have 
repeated the calculations for infinite systems with the same given parameter set. 
In Fig.l, the isotherms of symmetric nuclear matter (top panel) and asymmetric 
nuclear matter (bottom panel) with X =0.2 are displayed. The asymmetry pa- 
rameter is defined as X = (p n — p P )/ (Pn + P P )- The abscissa refers to po/p (= V/Vo) 
where po is the saturation density of normal nuclear matter given by 0.15 fm~ 3 . 
The isotherms are calculated with the NL3 parameter set and are shown for a few 
temperatures, at and around the critical temperature T c . For symmetric matter, 
T c is 14.2 MeV whereas for the asymmetric matter considered, it is 13.6 MeV. The 
isotherms resemble closely those obtained for the Van der Waals systems and are 
not quantitatively very different from those found in the nonrelativistic approach 
[§] . The dashed and the dotted lines in the figure are the liquid-gas coexistence line 
and the spinodal line, respectively. With the other parameter sets, the results are 
very similar; however, the critical temperature is found to increase and the critical 



volume becomes smaller for the parameter set which gives larger incompressibility. 
The critical temperature and the critical volume for the different parameter sets 
are given in Table I. 

For the symmetric infinite system, the liquid-gas coexistence line is obtained 
from the Maxwell construction as it is effectively a one-component system. Here 
the pressure and the chemical potential in the two phases are the same at a fixed 
temperature (Gibbs criteria) throughout their coexistence for any value of A where 
A is the liquid volume fraction (the gas volume fraction is 1 — A). Then the neutron- 
proton asymmetry is zero in both phases for all values of A. The asymmetric 
nuclear matter is a two-component system; here for any A, not only that the Gibbs 
criteria are to be fulfilled for thermodynamical coexistence of the liquid and gas 
phases, the overall asymmetry (neutron-proton ratio) has to be conserved which 
introduces added complications. Now, the pressure, chemical potentials and the 
neutron-proton ratio in both phases are in general changing functions of the volume 
fraction 0. We have taken these factors into account to determine the liquid-gas 
coexistence region. The spinodals shown in Fig.l are the isothermal spinodals 
referring to mechanical instability. The diffusive spinodal JL7| for the asymmetric 
matter is not shown here. 

Finite nuclei (even symmetric ones) with Coulomb interaction switched on 
behave like two-component systems. However, the ideas expounded earlier for 
asymmetric nuclear matter for the construction of the liquid-gas coexistence lines 
can not be employed to the case of finite systems. The density in finite nuclei 
is not uniform unlike infinite nuclear matter which is homogeneous in either the 
liquid or the gas phase. For asymmetric infinite system, we have found that for 
the liquid-gas coexistence, in general the neutron-proton ratio in the gas phase 
is much larger than that in the liquid phase. There is thus a phase separation 
between the neutrons and protons. For the construction of the isotherms for finite 
nuclei, the system is enclosed in a finite volume; then the self-consistent solution 
of the density profiles does not allow any significant neutron-proton phase sepa- 
ration because of the strong attractive unlike pair interactions compared to the 



interaction among the like pairs. Indeed, from numerical calculations, we find that 
the N/Z ratio throughout the nuclear volume is nearly the same (not varying by 
more than 10%) at all temperatures beyond T ~ 3.0 MeV. For a two-component 
thermodynamic system, liquid-gas coexistence occurs along the Maxwell-line (con- 
stant pressure, equal areas in the unstable phases in the P — V diagram) if the 
ratio of the concentrations of the components are the same in the two phases as 
it is then effectively a one-component system. It is known that for small finite 
one-component systems, the pressure in the coexistence region may have a small 



negative slope |18j in the P — V plane, but for symmetric nuclei with Coulomb 
switched off, we find numerically that the conventional Maxwell construction is 
an excellent approximation as the differences in the chemical potentials on both 
ends of the Maxwell line are negligibly small. It is found that for the asymmetric 
finite nuclei under study (even with Coulomb on), the difference in the neutron 
or the proton chemical potentials at the ends of the Maxwell line is typically ~ 
0.2 MeV only which is around 30-40 times smaller compared to that for infinite 
nuclear matter with the same asymmetry. We therefore expect that for these finite 
systems, for the determination of the liquid-gas coexistence, conventional Maxwell 
construction may not be a poor approximation to which we have resorted to in the 
present calculations for the nuclei considered. 

The isotherms for the lightest nucleus 40 Ca and the heaviest nucleus 150 Sm 
that we consider are shown in figures 2 and 3, respectively for the parameter set 
NL3. The results with the other parameter sets are not displayed as they look 
very similar. The finite size effects and the Coulomb interactions between protons 
do not change the qualitative character of the isotherms, the only effects are the 
lowering of the critical temperature and raising of the critical volume to some 
extent. The coexistence lines and the isothermal spinodal lines are shown by the 
dashed and the dotted lines, respectively. The critical temperatures and the critical 
volumes for different parameter sets are shown in Table I. The finite size effects as 
well as the Coulomb interaction tend to reduce the critical temperature. To isolate 
the Coulomb effect, the critical parameters for l50 Sm are also displayed in Table 



I switching off the Coulomb interaction. It is seen that the Coulomb interaction 
lowers T c for 150 Sm by about 1 MeV. The nature of EOS shows that it is possible 
to have a liquid-gas phase transition in a finite nuclear system below the critical 
temperature T c if it is prepared suitably in thermodynamic equilibrium. 

To make a quantitative comparison of the results in the RTF theory with those 
obtained in the non-relativistic Thomas- Fermi (NRTF) framework ||, the NRTF 
results for the critical temperatures and the critical volumes for nuclear matter 
and for the nucleus 150 Sm are also given in Table I. The NRTF calculations were 
performed with a modified Seyler-Blanchard (SBM) effective interaction which 
gives ph 240 MeV. This lies in between those obtained with NL1 and NL3 
parameter sets. However, it is seen that the critical parameters for the SBM 
calculations are in between NL3 and NLSH for infinite nuclear matter and close 
to those obtained with the NL3 parameter set for finite systems. The nature of 
the EOS and the critical parameters are controlled by the single-particle potential 
and the effective mass. To compare these quantities in the RTF and in the NRTF 
models, we display in Fig. 4 the nucleon single-particle potential as a function 
of the nuclear density (scaled with the normal nuclear density po) f° r symmetric 
nuclear matter as given in the two models. The single particle potential in the 
RTF is taken to be V + 5 as given by Eqs.(7) and (8). The corresponding effective 
masses are shown in Fig. 5. It is seen that the SBM single-particle potential is very 
close to that obtained with the NL1 parameter set at lower densities which evolves 
towards that generated with the NL3 parameter set with increasing density. The 
effective mass in the NRTF model is only a few percent lower at low densities and 
near the normal nuclear matter density it becomes somewhat higher compared to 
those obtained in the RTF model. The single particle potential for protons as 
a function of distance from the centre of the nucleus ( 150 Sm) in the two models 
corresponding to the ground state and at a temperature T = 8 MeV are shown in 
the top and the bottom panels of Fig. 6, respectively. The neutron single-particle 
potentials (not shown) have very similar behavior. It is seen that the potentials 
obtained in the NRTF model and those obtained with the different parameter sets 



in the RTF models are quite close. The effective mass for finite nuclei in the two 
models is consistent with that shown in Fig. 5, i.e., the NRTF model yields effective 
mass which is a little higher at the center and lower at the surface compared to 
the RTF ones. The above results clearly demonstrate the closeness of the single- 
particle potential and the effective mass in the non-relativistic and in the relativistic 
framework. So it is natural to expect that the EOS and the related thermodynamic 
properties would not be very different. 

3.2 The caloric curve 

We have remarked before that the density and hence the observables depend on 
the volume in which the nucleus at finite temperature is confined. The calcula- 
tion of the excitation energy as a function of temperature (the caloric curve) is 
thus volume dependent. In the experimental situation with energetic heavy ion 
collision, it is generally assumed that the hot nuclear system prepared after the 
collision expands substantially beyond its normal size (~ 4 — 8 times Vo) and then 
undergoes fragmentation due to density instabilities. Guided by the practice that 
many theoretical calculations for heavy ion collisions are done by imposing that 
thermalization occurs in a freeze-out volume, we fix a volume and then find the 
caloric curve. The freeze-out volume V is determined by r max occuring in Eqs. ( p6| ) 
- ( p3|) as V/Vq = (l/A)(r max /r ) 3 . Here r is the radius parameter corresponding 
to the normal nuclear volume Vo; it is taken to be 1.2 fm. 

In order to see the signature of liquid-gas phase transition, it is evident that 
the freeze-out volume is to be chosen beyond the critical volume V c . It is seen from 
Table I that for the nuclei considered, V c ~ 5 — 6 times Vo- We have fixed the 
freeze-out volume as 8Vo for our calculations. In Fig. 7, the caloric curves for the 
nuclei A0 Ca, 109 Ag and 150 Sm are displayed with the three different parameter sets. 
The caloric curves are seen to be nearly independent of the parameter sets. At 
lower temperatures , the excitation energy per particle e* increases quadratically 
with temperature similar to that in a Fermi gas; for T between 5-10 MeV, the 
caloric curve exhibits a shoulder and beyond that with a kink e* rises linearly with 



temperature as observed for a classical gas. In the bottom panel of Fig. 7 the 
data for experimental caloric curves for the ALADIN and the EOS collaboration 
are also shown. It is noted that the calculated caloric curves show shoulders at 
temperatures significantly higher than those obtained from experiments. This 
may be attributed to (i) the neglect of fluctuations in the theory which is expected 
to play important role near the transition temperature and (ii) the neglect of 
collective flow which lowers the transition temperature appreciably as seen in the 
non-relativistic calculations [§]. The specific heat C v defined by Eq. fl34f ) for the 
three parameter sets are shown in Figs. 8 and 9 for the three systems studied. 
Except for the top pannel in Fig. 8 , the calculations reported are done at the 
freeze-out volume 8Vo- The heat capacity shows a sharp peak signalling the liquid- 
gas phase transition, the peaks occurring at those temperatures where the caloric 
curve exhibits a kink. It is found that the harder the EOS the larger the transition 
temperature (T p ) and the sharper the peak. In the top panel of Fig. 8 , the specific 
heat for 40 C a with a freeze-out volume 4Vo is displayed. Instead of a sharp peak, 
it shows a broad bump at T ~ 10 MeV. Here the freeze-out volume is less than the 
critical volume, we do not associate this bump with a liquid-gas phase transition; 
it possibly signals a precursor to the transition. 

The evolution of the density distributions for the three systems around the 
phase transition temperature T p are displayed in Figs. 10 and 11. The calculations 
are done in a freeze-out volume 8Vo- In Fig. 10, the density distribution for 40 Ca 
with the parameter sets NL3 (top panel) and NLSH (bottom panel) are shown. At 
low temperatures, the density is more like a Woods-Saxon profile; with increasing 
temperature, the central density depletes and a long tail spreading to the boundary 
develops as is shown by the dashed line at T = 0.9T P . With further increase in 
temperature to T p and a little beyond, the density distribution tends to be uniform 
as is evident from the solid line (T = T p ) and the dotted line (T = 1.1T P ). It is 
found that the evolution of density with temperature is nearly independent of 
the parameter sets and therefore for the nuclei 109 Ag and 150 Sm calculations with 
only NL3 parameter sets are shown in Fig. 11. For these nuclei too the evolution of 



density is very similar to that in 40 C a. The rapid change in the density distribution 
towards a uniform one as the temperature approaches T p is a further indication of 
a liquid-gas phase transition in these finite systems. 

For symmetric nuclear matter, at a fixed temperature, the pressure remains 
constant over the whole coexistence region. The Gibbs free energy per particle g 
then shows a kink at the transition temperature when the pressure is held con- 
stant. Its derivative with respect to temperature, the entropy function, then shows 
a discontinuity there. For asymmetric nuclear matter, it has been shown in Ref. 
that the liquid-gas phase transition is second-order, the continuous transition be- 
coming more conspicuous with increasing asymmetry. There the pressure is not 
constant but shows a negative slope in the coexistence region for an isotherm in 
the P — V plane. Thus at constant pressure, the end points of the coexistence line 
are at different temperatures, then the kink in the Gibbs free energy g disappears 
and the entropy function becomes continuous with kinks at the end points of the 
transition. The liquid-gas phase transition thus occurs over a finite temperature 
interval. For finite nuclei, the exact calculation of the thermodynamic functions in 
the coexistence region is nontrivial and still not known, but as explained earlier in 
the beginning of this section, we expect the pressure in this region to remain close 
to a constant; it is then obvious that the entropy at constant pressure would show a 
discontinuity at a transition temperature as shown in Fig. 12 indicating a first-order 
phase transition. Precise statement about the order of the phase transition in finite 
nuclei can be made only if the exact calculations of the partition functions and 
hence the thermodynamic variables are rendered posssible. The present mean-field 
theory shows that the liquid-gas phase transition in finite nuclei has characteristics 
closer to a first-order transition. This is in consonance with the results obtained in 



the lattice gas model [fL9fl . No definite conclusion can however be reached from the 
experimental data; the GSI data Jl| indicate a first-order phase transition whereas 
the analyses of the EOS data |2T| show that the transition may be second-order. 



4 Conclusions 



The relativistic mean-field theory which has been very successful in describing the 
ground state properties of nuclear systems has been applied for the first time in 
evaluating the equation of state of finite nuclei in this paper. We have resorted to 
a relativistic self-consistent Thomas-Fermi theory with the non-linear a — u — p 
version of the Lagrangian. It is found that the results do not differ qualitatively 
from those obtained in a non-relativistic approach. This is due to the fact that the 
single particle potential and the effective mass which control the relevant observ- 
ables are very similar. The calculations have been done with three parameter sets 
with very different nuclear incompressibilities, still the EOS look nearly the same 
and the critical parameters are not too different. The critical parameters for finite 
nuclei are appreciably different from those of the infinite system. The specific heat 
C v calculated from the caloric curve shows peaked structures signalling liquid-gas 
phase transition in the nuclei studied. The near uniformity of the density distri- 
bution as the system approaches the transition temperature confirms this further. 
Analysis of the thermodynamical quantities indicates that this liquid-gas phase 
transition in the finite nuclei is more compatible with a first-order one. 
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Table 1: Critical temperature and critical volume for a few systems in the RTF 
model with parameter sets NL1, NL3 and NLSH; and in the NRTF model with 
the SBM interaction. 



Systems 


T c (MeV) 


Vc/Vq 


NL1 


NL3 


NLSH 


SBM 


NL1 


NL3 


NLSH 


SBM 


Sym. NM 


13.4 


14.2 


15.4 


14.5 


3.9 


3.3 


2.9 


2.8 


Asy. NM 
(X = 0.2) 


12.7 


13.6 


14.7 


14.1 


4.1 


3.4 


3.0 


2.9 


m Ca 


11.1 


11.6 


12.4 




6.7 


6.5 


5.8 




150 Sm 
(With Coulomb) 


11.4 


12.0 


13.0 


11.8 


5.4 


5.0 


4.4 


4.9 


150 Sm 

(No Coulomb) 


12.3 


12.9 


14.0 


12.5 


6.4 


6.0 


5.2 


6.2 



Figure Captions 

Fig. 1 The equation of state for symmetric nuclear matter (top panel) and of asym- 
metric nuclear matter (bottom panel) with NL3 parameter set. The temper- 
atures (in MeV) for the isotherms are as marked in the figure. The dotted 
lines are the spinodals and the dashed lines are the coexistence curves. 

Fig. 2 The equation of state for the nucleus 40 Ca for the NL3 parameter set. The 
different notations used have the same meaning as in Fig. 1. 

Fig. 3 Same as in Fig. 2 but for the nucleus 150 SVn. 

Fig. 4 Plot for the relativistic and non-relativistic mean field potentials as a function 
of density for symmetric nuclear matter at zero-temperature. 

Fig. 5 Density dependence of nucleon effective mass for symmetric nuclear matter 
at zero temperature calculated within the relativistic and non-relativistic 
framework. 

Fig. 6 Relativistic and non-relativistic mean field potentials for protons at tempera- 
tures T = MeV (top panel) and T = 8 MeV (bottom panel) for the nucleus 
150 Sm. 

Fig. 7 Caloric curves for the systems 40 Ca, 109 Ag and 150 Sm. Open circles and filled 
squares in the bottom panel represent the experimental data for ALADIN 
and EOS collaborations. 

Fig. 8 The specific heat for 40 C a at freeze-out volume equal to 4Vo (upper panel) 
and 8Vb (lower panel) with different parameter sets as labelled. 

Fig. 9 The specific heat for the nucleus 109 Ag (upper panel) and 150 Sm (lower panel) 
at the freeze-out volume 8Vo with different parameter sets. 

Fig. 10 Density distributions around transition temperatures T p = 10.1 and 10.3 
MeV at a freeze-out volume 8V for m Ca with parameter sets NL3 (upper 
panel) and NLSH (lower panel), respectively. 



Fig. 11 Density distribution around the transition temperatures T p = 10.2 and 10.3 
MeV for 109 Ag (upper panel) and 150 Sm (lower panel), respectively with NL3 
parameter set at a freeze-out volume 8Vo- 

Fig. 12 The temperature evolution of entropy per particle for the system 150 Sm at a 
constant pressure (P =0.087 MeV/fm 3 ). The dashed line shows the discon- 
tinuity at the transition temperature. 
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